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Abstract 



It is shown how the use of the Kerr-Schild coordinate system can greatly simplify the formula- 
tion of the geodesic equation of the Schwarzschild solution. An application of this formulation 
to the numerical computation of the aspect of a non-rotating black hole is presented. The 
generalization to the case of the Kerr solution is presented too. 



PACS numbers: 0230H, 0470B 



1 Introduction. 



Since the publication of the Schwarzschild static solution solution , it has been well known that the 
geodesic equation in this space-time can be solved analyticaly. The purpose of this note is to show 
that, by the mean of the Eddington coordinate system, the equation of motion of test particles in the 
Schwarzschild space-time can be greatly simplified so as to provide a more efficient method of solution, 
using ideas suggested by analysis of the less simple case of the Kerr rotating black hole solution. 

The forms of the spherical Schwarzschild metric as expressed in term of the usual Schwarzschild 
coordinate system and in term of the less usual Eddington 1 1] coordinate system (which is interpretable 
as a limiting case of the Kerr-Schild § coordinate system for a rotating black hole) are presented in § 2. 
Taking advantage of the obvious constants of motion which appear in the Schwarzschild space-time, we 
write down the geodesic equation in term of the Eddington coordinate system as second order equations 
with respect to some affine parameter. 

These results will be useful for many problems for which geodesic motion is concerned. In particular, 
they can greatly simplify the numerical integration of the test particles trajectory. As an application, we 
present in § 3 the apparent shape of a non-rotating Black Hole surrounded by an accretion disc as seen 
by an observer who is travelling towards and ultimately through the event horizon. 

2 The Schwarzschild metric and the geodesic equation. 

Using the standard Schwarzschild cordinate system (t,r,9,(p) 7 the line element of the metric takes the 
form 

ds 2 = -(1 - — )dt 2 + \rrdr 2 + r 2 (d6 2 + sin 2 6dtp 2 ) (1) 
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in units such that c = G = 1, where M is the mass of the Black Hole. 

The first integrals of the equations of motion are well known to be expressible as 
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where a dot denotes differentiation with respect to some affine parameter, r say, and where the constants 
of motion a 2 , E and <E> are respectively the rest-mass, the energy and the angular momentum about the 
axis sin 6 = of the test-particle, K being the Carter's fourth constant of motion which reduces in 
this simple spherical case to the square of the total angular momentum of the orbiting particle. 

Introducing the Eddington coordinate system (T, x, y, z), i.e. the Kerr-Schild coordinate system for 
the limiting case of a non rotating black hole, where (x, y, z) are "cartesian-like" spatial coordinates and 
where T is a retarded time, which are related to the Schwarzschild cordinates by means of 
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and 



x = r sm v cos tp; y = r sin 6 sin ip; z — r cos ( 
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T = u — r ; du = dt 
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At this stage, it is to be noticed that the above coordinate transformation (first published by Eddington 
in 1924 solves the "Schwarzschild singularity" (r = 2M) problem in the sense of being regular accross 
the event horizon r = 2M . 

In terms of this system, it is straightforward to show that the geodesic equations reduce to 



the metric takes the (Kerr-Schild) form 

ds 2 = -dT 2 + dx 2 + dy 2 + dz 2 

where 



(xdx + ydy + zdz + rdT) 
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r = x + y + z 
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(7) 


x l = -3KM^ - 




(8) 



where x l := (x,y,z). 

Keeping in mind that \f\ = E when any test-particle reaches the event horizon, the above equations 
show clearly the irreversible nature of the Schwarzschild Black Hole. Moreover, unlike what is obtained 
in terms of the Schwarzschild cordinates system, for which one has to take care of the axis sin# = 0, of 
the sign of r and 6 and of the presence of the event horizon, such differential equations can be numericaly 
solved in a straightforward way in order to determine the path of any test-particle. 



3 Apparent shape of a non-rotating Black Hole. 

Accretion disc are currently supposed to play an important role in several astrophysical situations, espe- 
cially when high-energy phenomena are involved. Because of the astrophysical interest of such objects, 
several authors have computed the apparent shape of of black hole with thin accretion disc as seen from 
infinity ([0 |B[ |^|). The main purpose of these calculation was to simulate line profiles which are com- 
monly observed from astrophysical sources which are generally interpreted as emission from an accretion 
disc around a compact object which may be a black hole. Our present purpose is quite different in the 
sense that, as an application of the previous formulae, we will present the apparent shape of a black hole 
surrounded by a thin accretion disc as seen by an observer who is flying near the hole. 

We will assume that the disc is a stationnary Keplerian one, orbiting in the plane z = 0. Each particle 
of the disc then follows a circular orbit with angular velocity 
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Following Page and Thorne pi , the flux of radiation from the surface of the disc is given by 
where M is the accretion rate and where we have introduced 

»=\fi- < n > 

As shown by Ellis || , the observed bolometric flux F D is given by 

F ° = JT+W (12) 

where the redshift factor (1 + z) is given by the ratio of the two scalar products 

\ PctW \emission /1Q\ 

(1 + Z) — — , \Y6) 
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p, w and u being the four-velocities of the photon, the emitting particle and the observer respectively. 
This redshift factor consists of a gravitational part due to the gravitational field of the black hole (which 
is measurable only in the closed vicinity of the hole), a Doppler part due to the rotation of the disc (the 
dominant one) and a Doppler part due to the motion of the observer. 

It is now straightforward to compute how the black hole will look. Assuming that the observer takes 
a camera with him, the coordinates of each pixel of the plate of the camera allows to compute the 
components of the four-velocity of the photon that reaches that pixel (recall that the four-velocity of a 
photon has only two true degrees of freedom) . The four velocity of each photon reaching the eyes of the 
observer can be obtained in terms of an orthonormal parallel-propagated frame (A , Aj , A 2 , A 3 ) along 
the observer's trajectory as 

p^ = Aq + cos "dXi^ + sin •& cos tp\^ + sin $ sin y>Ag , (14) 

where I? and if are two angles describing the photographic plate in the rest frame of the observer. 

A direct numerical integration of equations ([?]) (||) in the special case [i 2 = for negative values 
of the affine parameter r gives the history of the photon and then the luminosity of the source of its 
origin. Finally, application of the correction factor (1 + z) 4 as given by equations ( |l2| ) and ( |T3| ) gives the 
brightness of this fiducial pixel. 



The figures represent eight simulated photographs obtained by this procedure. These "snapshots" have 
been computed at succesive steps during the flight of an observer who is on free fall along a parabolic 
orbit (/i 2 = 1 , E = 1, K = 12) of the Schwarzschild space-time (see figure 1). On the first snapshop, 
the observer is above the disc at a distant r = 1280M. The observer then crosses the disc (snapshot # 4, 
r = 39M) and, finally, goes into the black hole (snapshot #8, r — 0.7 M). During this trip, the observer 
directs his camera in the direction of his motion except on the last snapshot for which the observer directs 
his camera toward the exterior of the hole. 

The brightness of the pictures has been computed by means of formulae © and ^) while the 
coloration has been added artificially. The apparent position of distant stars has also been computed. 
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snapshot #6 



Figure 1 : schematical representation of the observer's path. 

This helps to show the gravitational lense effect which is more conspicuous on snapshot # 2 (r = 121 AT): 
some distant stars, which are more or less perfectly aligned with the observer-black hole system, appear 
as pieces of rings in the sky. The opacity of the disc has been extremely minimized in order to make 
more visible the second (and even third) image of the disc and the distant stars. Moreover, this allows 
to show both the apparence of the upper and down sides of the disc. 

The trajectories of the photons and the link between these trajectories and the apparent shape of 
the hole as seen by the observer are schematically given in figure 2. Examination of figure 2 helps 
to understand the apparent shape of the disc. Some of the photons emitted by the disc are strongly 
deflected by the gravitational field (not the curvature) before reaching the eyes of the observer. Hence, 
the observer is able to see the upper rear side of the disc which appears like surrouding the top of the 
black hole (primary image on figure 2). In the same way, the observer can see the down rear side of the 
disc which seems to surround the bottom of the black hole. Now, some photons whose impact parameter 
is very close to the capturing one make one (even two) turn(s) around the black hole before reaching the 
eyes of the observer. This leads to the formation of a second (third) image of the disc. 

The numerical scheme used in these calculation is of second order. The integration has been done 
where respect to the parameter A which is related to the affine parameter by 

<*A = % (15) 

and which is well suited to the form of equations (0-||) in the sense that it exploits the fact that, far from 
the hole, null geodesies are straight lines. The code is fully parallelized. Vectorisation is also possible but 
less efficient since the number and the kind of operations needed to compute the value of each pixel of 
the screen vary from a pixel to another one. 
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Figure 2 : schematical representation of the photon trajectories. 



4 Conclusion. 



We have shown that the use of the Eddington coordinates system in the Schwarzschild space-time facil- 
itates the description of the geodesical motion in that space-time. This formulation is well adapted to 
direct numerical integration because it avoid the usual troubles tied to the spherical-type coordinates 
system and the trouble tied to the pseudo-singularity r = 2M. 

Such a formulation can be extended to the non-spherical case of the Kerr black hole solution. One 
can show that, in the Kerr-Schild coordinate system, the equations of motion of a test particle of zero 
rest-mass can be written in the form: 

Y?— = -4or^£j (16) 
M A 

-(£ - 4r 2 ) sin cos -0 - ^a- 

-ar sin 9 sin i/jQ j4(£E -Q) + (4a 2 - S)^ 

S 3A = +AaS^x (17) 
M A 
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+ar sin cos ip-^ < 4(££ - Q) + (4a 2 



E 3 -^ = -Xcos^Sr 2 -a 2 cos 2 (9) (18) 
= E(^)^ (19) 



where we have introduced the quantities 



£ = r 2 + a 2 cos 2 , (20) 
£ = E(r 2 + a 2 ) 2 - a$ , (21) 
Q = £r + £ (22) 

These results have been used in the study of profiles and shifts of lines emitted from Keplerian accretion 
disc around in X-ray binaries containing a neutron star or a black hole H and to the study of microlensing 
effects in active galactic nuclei M . 
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